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Abstract 

Misperceptions about extreme dependencies between different financial assets have been an im- 
portant element of the recent financial crisis. This paper studies inhomogeneity in dependence 
structures using Markov switching regular vine copulas. These account for asymmetric depen- 
dencies and tail dependencies in high dimensional data. We develop methods for fast maximum 
likelihood as well as Bayesian inference. Our algorithms are validated in simulations and applied 
to financial data. We find that regime switches are present in the dependence structure of various 
data sets and show that regime switching models could provide tools for the accurate description 
of inhomogeneity during times of crisis. 
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1. Introduction 



It is a well-known fact for univariate financial time series, that variances are not constant over 



time and tend to cluster together, as it is described in Engle's seminal paper (Engle 1982). In the 



multivariate setting, however, investigating whether also the dependence structure between different 
assets is varying over time is far more challenging. 

In this paper, we propose to combine the very flexible class of regular vine (R-vine) copulas with a 
hidden Markov structure, which accounts for the stylized facts observed in the dependence structure 
of multivariate financial return series and to investigate whether different dependence periods are 
present in given data. It is observed for economic data that the dependence is asymmetric between 
negative and positive returns and that also tail dependencies exhibit asymmetries, i.e. differ between 



the upper and lower tail (Longin and Solnik (2001) , Ang and Chen (2002) |. While classical models 
for multivariate time series, which are based on multivariate normal or Student-t distributions, fail 
to explain these observations, R-vine copulas utilize the richness of the class of bivariate copulas to 



account for such types of dependence (Joe et al. 2010). These copula models, which are constructed 



based on a series of linked trees called R-vine, constitute a new class of multivariate dependence 



models which has been introduced by Bedford and Cooke (2001 2002) using and generalizing ideas 



of Joe (1996) Since R-vine distributions are built up hierarchically from bivariate copulas as 
building blocks, also inference can be performed exploiting that hierarchical structure as pioneered 
in the seminal paper of Aas et al. (2009)] Further, they allow to scale dependence modeling to 



large dimensions by using truncation techniques as introduced by Heinen and Valdesogo (2009b) 



and Valdesogo (2009) and considered more general in Brechmann et al. (2011)] The hierarchical 
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structure is exploited for very fast yet asymptotically efficient sequential estimation (Half 2010). 



Applications in dimensions as high as 52 (Brechmann and Czado 2011 1 or 100 (Heinen and Valdesogo 



2009a) have been considered using restricted subclasses of the general R-vine model we propose. 



Being exposed to the challenge of high dimensionality in many applications, research in the area of 
multivariate dependence modeling is focussed on the case of time-homogeneous dependence struc- 
tures. However, there are also promising approaches for allowing variations in the dependence 
structures over time. Hereby, one popular direction of research uses parametric dependence models 



in which the parameters constitute a function of time. For example, Hafner and Manner (2010) 



and Almeida and Czado (2011) model dependence parameters as a latent autoregressive stochastic 
process. Another popular direction combines multivariate dependence models with regime switching 



and in particular Markov switching (MS) models. Publications in this direction include Pelletier 



(2006)1 [Garcia and Tsafack (2011)[ [Okimoto (2008) | [Markwat et al. (2009)1 and in particular 



ChoUete et al. (2009) who first applied a restrictive vine distribution in the context of time vary- 



ing dependence models. A recent survey, comparing these different approaches with focus on the 



bivariate case, is given in Manner and Reznikova (2011) 



In this paper, we will pick up the second approach and extend the applicability of regime switching 
dependence models to high dimensions using a general R-vine model, extending the initial approach 



of ChoUete et al. (2009) Since high dimensionality and the presence of latent state variables make 
parameter estimation challenging and much of econometric literature on hidden Markov and MS 
models only marginally considers inference methods, our main concern will be to fill this gap. 
In this respect, we give a thorough introduction to parameter inference for MS models, extending 
the existing procedures in several ways. While existing models for regime switching dependence 
structures lack the possibility to be extended to high dimensions, we demonstrate how parameter 
inference for the MS R-vine model can be facilitated in almost arbitrary dimension by using a 
fast approximative Expectation - Maximization (EM) procedure in a Maximum Likelihood (ML) 
framework. The developed algorithm can also be used to obtain proper starting values to per- 
form Bayesian parameter inference using MCMC techniques. We extend the Bayesian estimation 



procedure which has been developed by Min and Czado (2010) for Student-t copulas on drawable 
(D-)vines to the case of general R-vine distributions with arbitrary bivariate copulas and include 



inference about the underlying MS model as it has been considered by Kim and Nelson (1998) 



This Bayesian inference method enables us to assess the uncertainty in parameter estimates by con- 
sidering Bayesian credible intervals (CIs). This has not been possible for MS dependence models, 
because of the computational burden necessary to achieve bootstrapped confidence intervals. 
In order to demonstrate the applicability and performance of our procedures for parameter esti- 
mation, we perform a simulation study and investigate several applications to empirical data. In 
this context, we will also show how model selection can be performed for time- varying dependence 
structures conducting a rolling window analysis and compare different models using the Bayesian 
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deviance information criterion (DIG, Spiegelhalter et al. (2002)). 



The remainder of this paper is structured as follows: Section [2] introduces the Markov switching 



regular vine copula model by first introducing R-vine distributions in Section 2.1 and then combining 
them with an underlying Markov structure in Section 2.2 Section [3] then focusses on the problem of 
parameter estimation, discussing the stepwise procedure in Section |3.1| and Bayesian estimation in 
Section [3. 2 1 In Section [4] we present the results of our simulation study before turning to empirical 
applications in Section [5] We consider three data sets: nine exchange rates in Section 5.1 five 
Eurozone country stock indices in Section [5. 2[ and the returns of ten selected stocks in the German 
stock index DAX in Section 5.3 Section [6] concludes our paper and gives an outlook to directions 
of future research. 



2. The Markov switching regular vine copula model 
2.1. Regular vine distributions 

Regular vines as a graph theoretic tool for the construction of multivariate distributions have been 



introduced by Bedford and Gooke (2001 2002). The R-vine itself is specified as a sequence of 



linked trees. It can be used to construct a multivariate distribution model by assigning to each 
edge a copula which corresponds to a bivariate conditional distribution, determined by the vine. 
For this construction to be possible, the regular vine V on d variables, which consists of a sequence 
of connected trees Ti, . . . ,T(i-i, with nodes Ni and edges Ei, 1 < i < d — 1, needs to satisfy the 



following properties (Bedford and Gooke (2001) |: 



1 . Ti is a tree with nodes A''i = {1, . . . ,d} and edges Ei . 

2. For z > 2, Tj is a tree with nodes Ni = Ei^i and edges Ei. 

3. If two nodes in Tj+i are joined by an edge, the corresponding edges in Ti must share a common 
node, {proximity condition) 

There are two popular subclasses of R-vines which differ in the number of edges per node on each 
level. We call an R-vine 

• Canonical vine (G-vine) if in each Tree there is one node which has edges with all d — i 
other nodes. 

• Drawable vine (D-vine) if each node has edges with at most two other nodes. 

A five-dimensional R-vine is shown in Figure [T] The notation we employ throughout our paper 



follows Gzado (2010) To build up a statistical model using the R-vine, we associate to each edge 



j{e),k(e)\D(e) in Ei, for l<i<d — 1, a bivariate copula density Cj^g) fc(e)|£)(e)) i-^- the copula 
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Figure 1 An R-vine tree sequence in five dimensions witli edge indices. 



associated to an edge e is chosen to correspond to the copula of the conditional distribution indexed 
t>y k{e)\D{e). We call j(e) and k{e) the conditioned set while D{e) is the conditioning set. Let 
^D{e) denote the subvector of a vector X determined by the set of indices D{e) and X_j the subvetor 
of X where the j*^ entry is deleted. For the definition of the regular vine distribution we associate 
the bivariate copulas of the vine with bivariate conditional distributions. Let X = {Xi, . . . , X^i) be 
a random vector with marginal densities /i , . . . , /^^ respectively, and copula densities corresponding 
to the conditional distributions of -'^j(e) and X;t(e) given X^^ equal to Cj(e),fc(e)|De) then we call its 
distribution a regular vine distribution. In this case, the density is uniquely determined and given 
by 

d d-1 

fl,...Axi, ■ ■ ■ ,Xd) =YlMxi) ■ Yl n Ci(e),fc(e)pe(^(a^i(e)|xDj,(F(Xfc(e)|xz)J) (1) 
i=l 1=1 e£Ei 



as shown by [Bedford and Cooke (2001) If the marginal densities are uniform on [0, 1], we call the 
distribution in ([T]) an R-vine copula. Given an R-vine V, a set of corresponding parametric bivariate 
copulas B and their parameter vector 0, we denote the R-vine copula density by c(.|V, B, 0). 
While also other iterative decompositions of a multivariate density into bivariate copulas and 
marginal densities are possible, R-vine distributions have the particularly appealing feature that 
the values for F(xj(e)|x-De) and F(xjrj(g)|x£)^) appearing in Equation (HI can be derived recursively 
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without high dimensional integrations. To see how the tree structure is exploited for this, we re- 
fer the reader to Appendix [A] where it will become clear in the context of step-wise parameter 
estimation. 

There are many possible R-vine structures in d dimensions and the indices in ([T]) do have a special 



order only in subclasses like D-vines or C-vines, as considered by Aas et al. (2009) For this reason, 
we apply the approach of Morales-Napoles (2008) Difimann (2010) and|Di£mann et al. (2011) to 



store the indices in a d x d lower triangular matrix M = {mij\i > j), where each row corresponds to 
a tree in the vine. In particular, given indices mjj, rn-j^j and (rrifc+i^j, . . . , rridi) in the matrix, there is 
an edge in the R-vine with {j(e), A;(e)} = {ma^mki} and D[e) = {mfc+i^j, . . . , nidi}. For example, 
the matrix corresponding to the R-vine in Figure [T] is 



/5 
2 
4 
3 



4/ 



The types of bivariate copulas in B as well as their parameters can conveniently be stored in 



matrices related to M, for details we refer to Difimann (2010) 



2.2. Markov switching models 

In this section, we introduce general MS models, as they have been established in statistics and 



econometrics by Hamilton (1989) focussing on the special case of multivariate dependence modeling. 



MS models constitute a special class of regime switching models, in which the process describing 
regime switches has a Markov structure. In essence, they assume that a hidden underlying process, 
which can be understood e.g. as the "state" of the world or the economy in financial applications, 
influences the development of a time series. 

In particular, we will assume that the underlying process affects the dependence structure of a 
multivariate time series. The different possible dependencies can then be described by different R- 
vine copulas for each regime /c, of which we assume the densities c(.| (V, B, 0)^) to be given. At each 
point in time, the present regime will determine the copula of the multivariate time series. For this, 
let {St)t=i,2,... be a homogeneous discrete time Markov chain with states {1, . . . For simplicity, 
we assume it to be of first order such that it can be completely characterized by its transition 
matrix P with elements P{St = i\St-i = j) = Pij- In applications where certain regressor variables 
influencing the development of the time series are known, our model could easily be extended to 
the case where the probabilities in this matrix change over time depending on these variables (c.f. 



Filardo and Gordon (1998) ). 
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Given the Markov chain {St)t=i,2,..., a MS R-vine copula model for a multivariate time series (uf = 
{uit, . . . , Udt),t = 1,2,...) where ua G [0, 1] Vi, t can be fully characterized by specifying conditional 
densities 

P 

c{ut\{V, B, St) = hst=k} ■ c(ut|(V, B, e)k). 

k=l 

The complete MS R-vine model is thus specified in terms of p R-vine copula specifications and the 
matrix P which contains the parameters of the underlying Markov chain. Considering inference in 
this context, we will always assume the R-vine structures Vk and corresponding sets of copulas B^, 
A; = 1, ... ,p, to be given and thus suppress them in the following notation. The MS R-vine copula 
is then completely described by its parameters 

where the subscript "cop" stands for copula parameters and "MC" for parameters making up the 
transition matrix P. Note that, although until now this is a pure copula model where no marginal 
time series structure is included, this introduces serial dependence. Given previous realizations it 
will be more or less likely that the hidden variable St assumes a specific state. The individual 
marginal time series {ui^t)t=i,2,... however are i.i.d. uniform for i = 1, . . . ,d. 



3. Inference for Markov switching models 



The first problem in developing inference methods for Markov Switching models is that we are faced 
with unobserved latent variables. In order to derive an expression for the full likelihood of a time 
series of observations ut = (ui, . . . , u^), let us consider a decomposition of their joint density into 
conditional densities: 

T 

f{nT\e) = f{u,\e).l[f{nt\nt-i,e) 



t=2 



Y,f{^i\Si = k,ek)P{Si = i\e 



MC) 



.k=l 



n 

t=2 



f{^t\St = k, Ok) ■ P{St = k\ut-i,e) 



.k=l 



where Uf := (ui, . . . ,Ut). The unconditional probabilities P{Si = i) in this expression are known 
from the stationary distribution of the Markov chain, which we assume to exist. To obtain the 
probabilities 

{nt\t-m)k-=PiSt = k\vit.ue) 



we can apply the filter of Hamilton (1989) Assuming Qit-i\t-i to be given, we calculate 



nt\t-i{e) = p ■ nt,i\t-ii0) and 

^t\t-i{0)Q{f{ut\St = k,ut-i,ek))k=i,...,p 



ELi fi^t\St = k, ut-i, e, 
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and obtain all probabilities, which are required to evaluate the density, recursively. The operator 
denotes componentwise multiplication of two vectors. 

Because of the latent state variables and the resulting dependence between parameters, direct 
maximization of the likelihood for given (Vfc,Bfc), k = 1, . . . ,p. is analytically not possible and 
numerically difficult. In the following, we discuss a frequentist and a Bayesian approach to make 
inference for this kind of model tractable. 

3.1. EM algorithm for MS models 



Hamilton (1990) proposed to overcome the problems in maximum likelihood estimation for an MS 
model by using an Expectation - Maximization (EM) type algorithm. This algorithm iteratively 
determines parameter estimates 6\ I = 1,2, . . . , which, under several regularity conditions, converge 
to the ML estimate for / — )• oo. In particular, it iterates the following steps: 

1. Expectation step: Obtain the conditional probabilities of the latent states = (5*1, . . . , St) 
given the current parameter vector 6\ i.e. P{St = St|uT,^')- 

2. Maximization step: Maximize a pseudo likelihood for 6^^^ , where the probability of being in 
a latent state St = St is replaced by the probability from step 1. 

More precisely, let us define the expected pseudo log likelihood function for 6^~^^, given the obser- 
vations viT and the current parameter estimate 6\ as 

Q(0'+i; UT, e') := / log (/(uy, St|0'+^)) /(ur, StW), (2) 

where Jg^ is short notation for summation over possible values of St, i.e. 

/st 5(St) = E"i=i • • • YTst=i diSi = si,...,St = St) for an arbitrary function g of St- 

With this specification, step 2 can be rewritten as 

2. Maximize Q{e^^^;UT,0^) with respect to 0'+^ 

We will now investigate how these steps can be performed in the context of the MS R-vine copula 
model. The notation in this section is chosen to follow the notation of Kim and Nelson (1999), for 
a justification of the algorithm we refer to the original work of Hamilton (1990). 
For simplicity, let us further assume that the Markov chain St has only two states, and denote 
a := Pii, b := P22, then 6mc = («> ^) and ^cop = (^1,^2)- Using the Markov property of St, we 
can decompose the joint density of (ut, St) as 

T T 

log (/(uT, St|0)) = log ifi^t\St, e,op)) + log {P{St\St-i, Omc)) ■ 
t=l t=l 
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With this, we can simphfy Expression n2K using 0^ = {d^copi^^Mc) 



log(/(uT,ST|0'+'))/(uT,ST|0') 



oc / log ( /(UT, St|0'+^) ) P(St|ut, 0') (with respect to 6') 



T 



J] log + J^log 

t=i t=i 

I log (/(ut|5i, e'+D) ■ p{St\ut, e') 



P(St|ut,0' 



+ 



St 



St 
T 



J^log P(5i|5i-l,0^+^, 



■P{ST\viT,e^) = A + B, 



where 



T 2 



T 2 



t=l St=l 



^ ^= E E (/(ut|5i,0i+i)) • p(St|uT,0') = 5^ E (/K|5t,0^+^)) • {^t\T{e'))s, 



t=l St=l 



and 



B :-- 



J^log (P(5i|5<_i,<+^, 



P{St\ut,0^)- 



Hereby, the probability (^ft^xi^'')) = P{St = St\uT,(^''), to which we will refer as the "smoothed" 
probability of being in state st at time t, can be determined from the output of the Hamilton filter 
by applying the following backward iterations, called Kim's smoothing algorithm. 



St 



where the division is to be understood componentwise. 

The second term B can be simplified further using the Markov property: 

T 



B = i j;iog(p(5i|5t_i,<+^)) 

-'^T \_t=l 



P(Sr|ur,0' 



T 1 



E (^('^* = = ^'mc)) nst = h St-i = i|uT, e') = 

t=l i,j=0 
T 

Y [{log{a)PiSt = 0, St-i = 0\UT, e') + log(l - a)P{St = 1, St^i = 0|ut, 9')+ 



t=i 



+ iogib)p{St = 1, St-i = i\uT, e') + iog(i - h)p[St = 0, St-i = i|uT, e% . 
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This has now to be maxhnized with respect to a and b. Taking the derivative with respect to a and 
b respectively, we obtain the step / + 1 estimates for a and b: 

j+i _T.J=iP{St = o,St.i = o\uT,e') 



a 



El=iP{St-i = o\nT,e') 
^ ELiP{St = i,St-i = i\nT,e^) 

In contrast to the model originally considered by Hamilton where all maximization steps could be 
performed analytically, this is not possible for the maximization of A with respect to the copula 
parameters in our case. This means that, while the transition probabilities can be obtained directly, 
the second part of the maximization step has to be performed using numerical optimization methods. 
Since a d-dimensional R-vine copula specification, in which each pair copula has k parameters, 
contains d{d — l)/2 ■ k parameters, this is computationally still very challenging. To circumvent 
this problem, we can exchange the joint maximization with the stepwise procedure of Algorithm [2] 
(Appendix |A|. 

We call this the stepwise EM- Algorithm. Since treewise estimation of copula parameters is asymp- 
totically consistent, this constitutes a close approximation to the "proper" EM-Algorithm. While 



there are theoretical results on the convergence of the EM-Algorithm (Wu 1983), we loose these 
properties with our approximation. All limit theorems however do rely on proper maximization at 
each step of the algorithm. This is almost impossible to guarantee in our case where we are faced 
with high dimensional optimization problems and have to rely on numerical techniques. Therefore, 
also an implementation of the "proper" EM-Algorithm has to be considered an approximation, 
which, apart from the asymptotic consistency of the step-wise estimator further justifies the use 
of a stepwise procedure. While all existing models for time varying dependence structures in high 
dimensions suffer from the computational burden for numerical estimation, we do only need to max- 
imize the likelihoods of bivariate copulas in this tree-wise procedure, which reduces computation 
time and avoids the curse of dimensionality. 

The estimate obtained by iterating the EM steps until convergence will be denoted by 

2^EM fp,EM ,p.EM -BM 
9 = [^cop =(^1 , - -,9k l^^MC 

3.2. Gibbs sampling for Markov switching models 

Having derived an approximative ML procedure for our MS models, we will now consider Bayesian 
estimation methods, which will enable us to quantify the uncertainty in parameter estimates. In 
particular, CIs and posterior standard deviations are determined naturally while the uncertainty in 
ML parameter estimates is very hard to assess in this context. Bayesian estimation for MS models 
has originally been considered by Albert and Chib (1993)] Building on their ideas, the Gibbs 



sampler which we develop in this section consists of updates for the copula parameters, the Markov 
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chain parameters and the latent state vector, respectively. To derive how these update steps are to 
be performed, let us first neglect the latent Markov part of the model specification and focus on the 
parameters 6 cop of R-vine copulas. 

Update of copula parameters 

In order to complete the model specification in a Bayesian framework, we first have to specify 



we assume non 



prior distribution for each component of Ocop- Following Min and Czado (2010) 
informative priors for all copula parameters in the model. For bivariate copula families considered 
where the parameter range is not compact, we restrict its support to some finite interval to avoid 
numerical instabilities for very small or large parameter values. If for all bivariate copulas there 
is a one-to-one correspondence between parameter values and Kendall's r given in closed form. 



the approach of Hofmann and Czado (2010) with uniform priors for r can be considered as an 



alternative. Furthermore, we can use a uniform prior for the correlation matrix of the model if all 



bivariate building blocks are Gaussian or Student copulas, c.f. Lewandowski et al. (2009) 
Since the posterior distribution is analytically not tractable and in particular expressions for the 
conditional posterior distributions are not available, we utilize the Metropolis-Hastings (MH) al- 
gorithm. As for the prior distribution, there are several possible choices for proposal distributions 



in this case. Min and Czado (2010) use a modification of standard random walk proposals where 
the normal distribution is truncated to the support of parameters, proposal variances are tuned to 
achieve suitable acceptance rates. This leads to poor acceptance rates in some cases with strong 
dependencies and to high autocorrelations in general. To overcome these problems, we consider a 
two point mixture of a random walk proposal with an independent normal distribution at the mode 
for each parameter. The modes of the individual parameters as well as their standard errors are 
approximated using the bivariate hessian from the stepwise estimation procedure for R-vines and 
both distributions are assigned a weight of 0.5. The methods of using independence proposals cen- 



tered around the mode has been proposed by Gamerman and Lopes (2006) and it has been applied 



in a context similar to ours by Czado et al. (2010) While there are parameter constellations where 



pure random walk proposals are more favorable than independence proposals and vice versa, simu- 
lation studies showed that the chosen mixture distribution works well for all settings. Given these 
considerations Step 3 of the Gibbs sampler can be performed using MH within Gibbs by sampling 
di for i = 1, . . . ,p from 

/(^j|u{ie{i,...T}|5t=i})- 
Update of Markov chain parameters 

For the second step of the Gibbs sampler, we will assume independent Dirichlet distributions as prior 
distributions for the columns of the transition matrix P, i.e. {Pi,j)i=i,...,n ~ Dirichlet((ajj)j=i^...^„) . 



11 



The likelihood function for the entries of P given a realization of the state vector St has the form 
of a multinomial distribution 

liP\ST) = flflpl'\ 
j=l i=l 

where riij denotes the number of transitions from state j to state i in St- Since the Dirichlet and the 



multinomial distribution are conjugate distributions (see Kotz et al. (2000)), also the conditional 
posterior distributions are Dirichlet distributions with parameters q,p°**<^"°'' = q^.^ _|_ j^jj-. From 
these we can sample directly. 

Update of the latent state vector 



For the first step, we follow the approach by Kim and Nelson (1998) who assume independent non 



informative priors for the latent states and draw St as a block. To do so, we decompose 

T-i r-1 
P(St|ut) = P(5t|ut) • n ^(^^l-^i+i'- • -^r'"^) = ^("^^l"^) ■ n ^(-^^l^^+i'"^)' 

1=1 i=l 

which allows us to generate St from P{St\ut) and St for t S {T — 1, . . . , 1} from 

P{St\ut,St+i) « P{St+i\St)P{St\ut), 

where P{St\ut) = Cl^t{6) can again be determined using the Hamilton filter. This corresponds to 
sampling St from /(St|ut, ^cop; ^A/c)- More formally, the three steps of the Gibbs sampler from 
above can now be reexpressed as follows: 

Algorithm 1 MCMC sampling for MS R-vine copulas. 
Sample St from /(St|ut,^) using Gibbs sampling. 
Sample P from /(P|St) using Gibbs sampling. 

Sample (0j|u|^g|]^ T}|5t=i})j fo^' i = ^, ■ ■ ■ ,n, using MH within Gibbs sampling. 

Iterating through these steps will yield 

^Qr,MCMC S^*'^'^^'^) = ^^(^^'.A/CAfC Qr,MCMC-^ Qr,MCMC^ St^^*^^*^) 

for r = 1, . . . ,R, where R is the number of sampled realizations from the posterior distribution. 

4. Simulation study 

Having derived all necessary components of our posterior sampling algorithm, this section sums up 
the results of a simulation study which has been performed in order to demonstrate the ability of 
the developed Bayesian inference procedure to capture the true model in simulated data. 
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We consider two regimes in four dimensions: A multivariate Gaussian copula, modeled as a D-vine 
copula where all bivariate copulas are Gaussian, and a C-vine copula, where all bivariate copulas are 
Gumbel or rotated Gumbel copulas; these regimes are kept fixed. We further consider 6 scenarios 
with different model parameters, summarized in Table [T] in terms of the corresponding values for 
conditional Kendall's r. In all scenarios we set the Markov parameters to a = 0.95 and b = 0.9, the 
corresponding prior distributions are chosen to be non informative. 



coverage probability 

conditional Kendall's r 90% CI 95% CI 







Gumbel regime 






Gauss regime 




symm. 


HPD 


symm. 


HPD 






T43|21 0.4 






T41|23 = 0.4 












Scenario 1 




r42|i = 0.6, T32|i = 0.6 






T42|3 = 0.6, r3i|2 = 0.6 




91.7% 


91.7% 


94.2% 


94.2% 




T41 


= 0.8, rai - 0.8, T21 - 


0.8 


T43 


= 0.8, T32 = 0.8, T21 = 


0.8 














■^43121 = 0.4 1 






T41123 = 0.1 












Scenario 2 




T4211 = 0.6, T3211 = 0.6 






T42|3 = 0.2, T3112 = 0.2 




89.2% 


89.2% 


91.7% 


91.7% 




T41 


= 0.8, r3i = 0.8, T21 = 


0.8 


7-43 


= 0.3, r32 = 0.3, T21 = 


0.3 














■''43121 = 0.1 






T41I23 = 0.4 












Scenario 3 




^4211 = 0.2, T32I1 = 0.2 






T42|3 = 0.6, r3i|2 = 0.6 




85% 


84% 


91.7% 


93.4% 






= 0.3, T31 = 0.3, T21 = 


0.3 


T43 


= 0.8, T32 = 0.8, T21 = 


0.8 














T43121 = 0.1 






T41123 = 0.1 












Scenario 4 




r42|i = 0.2, T3211 = 0.2 






T42|3 = 0.2, r3i|2 = 0.2 




75% 


75% 


81.7% 


92.5% 




T41 


= 0.3, r3i = 0.3, T21 = 


0.3 


T43 


= 0.3, r32 = 0.3, T21 = 


0.3 














^43121 = 0.3 






^41123 = 0.3 












Scenario 5 




T42|i = 0.5, T3211 = 0.3 






T42|3 = 0.5, r3i|2 = 0.3 




85.0% 


85.0% 


92.8% 


92.8% 




T41 


= 0.7, T31 = 0.5 T21 = 


0.3 


T43 


= 0.7, T32 = 0.5, T21 = 


0.3 











Table 1 Simulation scenarios investigated and empirical coverage probabilities based on 120 data 
sets from each scenario. 

From each scenario, we simulate a time series with 800 four dimensional observations. As an input 
for the MCMC algorithm, we need the R-vine tree structure of the model, the corresponding set 
of bivariate copula families, and starting values for the parameters. Keeping the (true) structure 
and copula families we used for simulations, we obtain a posterior estimate for the parameters as 
follows: 

1. Starting values for the EM algorithm: Fit the copula for each regime to the whole data set using 
the stepwise estimation procedure, and cluster the observations according to their likelihood 
values. Refit the copula to the 400 observations which have the highest log likelihood. Set the 
MC parameters to a = 6 = 0.9, since persistent regimes are expected. 

2. Starting values for MCMC: Iterate the stepwise EM algorithm until convergence. 
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3. Obtain 1000 independent samples from the posterior distribution of the parameters. The 
Gibbs sampler is started at the values obtained in Step 3, a burn-in period is discarded and 
the chain is sub-sampled according to the effective sample size. 

From the obtained samples, we estimate 90% and 95% symmetric and highest posterior density 
(HPD) CIs for the copula parameters and check whether all true copula parameters lie within these 
intervals. The HPD intervals are calculated using the "coda" package for R, while a% - symmetric 
intervals are determined by the empirical ^% and 1— ^% quantiles. The procedure was repeated 120 
times for each scenario with results reported in Table[T] Relative bias and mean squared error (MSE) 
of two selected scenarios are displayed in Appendix [B] Since the parameters are not independent, 
the combination of CIs for single parameters does not necessarily yield a credible region for the same 
level. This, in addition to the fact that Bayesian CIs in a simulation study are hard to interpret in 
terms of their frequentist coverage, makes a statistical analysis of the results difficult. 
But clearly the number of exceptions for each scenario, except Scenario 4, lies within the range of 
what we would expect, namely about 90% (95%) frequentist coverage. Scenario 4 corresponds to 
low dependence in both regimes, the higher number of exceptions in this case is due to identification 
problems. We conclude, that with clearly distinguishable regimes the outlined procedure is able to 
identify the true model. 



5. Applications 

In this section, we apply the MS model described in Section [2] together with the estimation proce- 
dures of Section [3] to analyze three financial data sets. Since the focus of this paper is on modeling 
dependence structures of multivariate data, we apply a two step estimation approach as suggested 
by Joe and Xu (199"6)] In the first step, appropriate parametric models for the marginal time series 



are fitted separately and used to transform the standardized residuals to approximately uniform 
margins. To this transformed data, we apply our copula model in the second step. While joint 
estimation of marginal and copula models as in Hofmann and Czado (2010)| is more effective and 



allows to take the estimation error in the models for the marginals into consideration, it is also 
computationally more challenging. Given the size of our data sets we choose the two step proce- 
dures. Note that the marginal time series structure we impose does not hurt the assumptions for 
our copula model, since the standardized residuals form approximately an i.i.d. sample. 
This proceeding implies that all statements regarding the dependence properties of investigated data 
do not refer to the dependence structure among the univariate marginal time series themselves. In- 
stead, we describe the dependencies of their transformed residuals obtained in the two-step approach. 
While the marginal models account for time-varying conditional variances and autocorrelations in 
given data, the dependence properties between individual variables are uniquely determined by 
the dependence properties of their residuals. This means that, although our approach does not 
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yield direct conclusions for observed variables, qualitative results and interpretations can always be 
transferred from the level of residuals to the level of observations. 

Before Bayesian or frequentist parameter inference for the MS R-vine model can be conducted, 
appropriate R-vine structures and sets of bivariate copulas for each regime need to be selected in 



a preanalysis. To do so, we apply the heuristic model selection techniques as outlined in Difimann 



(2010) Difimann et al. (2011) and Brechmann et al. (2011) These are based on the following 



steps: 



1. For each pair of variables estimate the corresponding value of Kendall's r from the copula 
data set. 

2. Create a fully connected graph which consists of the marginal variables as nodes and where 
an edge is added between every pair of variables. 

3. Associate to each edge the absolute value of the corresponding Kendall's r as a weight. 

4. Determine the maximum spanning tree (MST), i.e. find a tree which maximizes the sum of 
edge weights using for example the algorithm of Prim (1957)| 



5. For each edge in the resulting tree fit a parametric bivariate copula from a catalogue of 
bivariate copula familis and estimate its parameters by ML. The specific family is then chosen 
using the Akaike information criterion (AIC). 

6. For each pair of edges which share a common node apply the probability integral transforma- 
tion of the conditional cdf given the common node based on the copula parameters estimated 
in step 4 to the corresponding copula data. This respects the proximity condition and provides 
pseudo observations for the next tree. 

7. Proceed with the pseudo observations as in steps 1 to 6 until all trees together with their 
copula types and parameters are determined. 

In all our applications we assume the presence of two regimes. Unless mentioned otherwise, the 
copula families we will consider in step 5 are the Gauss copula and the Gumbel copula. Since 
the Gumbel copula is not invariant with respect to rotations, we consider its standard form and 
rotations by 90°, 180° and 270°, respectively. For all models studied, we run the MCMC for around 
20000 internal iterations discarding the first 1000 as burn-in, and keep every fifth observation to 
reduce autocorrelations. For estimating quantiles of the posterior distribution, we further thin the 



output according to what Kass et al. (1998) call the "effective sample size" (c.f. Carlin and Louis 



(2009) |. For stability reasons, we use this as a proxy for thinning with respect to the autocorrelation 



function itself. After this, we end up with ca. 1000 approximately i.i.d. samples. 
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5.1. Exchange rates 

The first data set taken into consideration consists of 9 excliange rates against tlie US dollar, namely 
the Euro, British pound, Canadian dollar, Australian dollar, Brazilian real, Japanese yen, Chinese 
yuan, Swiss franc and Indian rupee. The observed time period is from July 22, 2005 to July 17, 2009, 
resulting in 1007 daily observations. The modeling of the one dimensional margins with appropriate 
ARMA-GARCH models and the transformation to copula data has been performed by |Czado et al. 



(2011) In total, we consider 6 models, which will be defined as we proceed, their defining trees and 



the allowed copula families are listed in Table |2} 



Regime 1 
(no crisis) 



Regime 2 Copulas 
(crisis) Regime 1 



Copulas 
Regime 2 



Model (1) 


Vii 






V12 




mixed 


7^ mixed 


Model (2a) 


V21 - 


Vi 




V22 


- V2 


N 


7^ SG 


Model (2a*) 


V21 = 


Vi 




V22 


= V2 


N 


7^ SG, N 


Model (2b) 


V21 = 


Vi 




V22 


- V2 


N 


+ G 


Model (2c) 


V21 = 


Vi 




V22 


= V2 


N 


^ Student-t 


Model (3) 


V31 = 


Vi 




V32 


= V3 


mixed 


7^ mixed 



Table 2 R-vine models considered for the exchange rate data. 



5.1.1. R-Vine with switching parameters 

As a first model for the exchange rate data set, we consider an R-vine with only MS parameters. 
To do so, we fit an R-vine with corresponding bivariate copulas to the data using the outlined 
procedure allowing for Gaussian (n), Gumbel (G), survival Gumbel (rotation by 180°, abb. as SG), 
and 90° / 270° rotated Gumbel (G90 / G270). As the estimated parameters for the bivariate copulas 
corresponding to the higher trees indicate conditional independencies, we truncate the R-vine copula 
after the second tree, i.e. we associate all edges on higher trees with independence copulas. The 
R-vine copula structure resulting from this procedure is tabled in Appendix [C] We call this model 
Model (1). 

Figure [2] shows the probability P{St = 2\ut,0 ) that the hidden state variable St indicates the 
presence of Regime 2 plotted against time. While Regime 1 is predominant until around February 
2007, Regime 2 becomes more important during the later times of the financial crisis. 
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o 

2005-07-22 



2006-05-09 



2007-02-26 



2007-12-06 



2008-09-22 



2009-07-09 



Figure 2 Estimated probabilities over time of being in state 2 for Model (1). (Solid: EM estimates 
smoothed by an MA(7) filter, Dotted: Bayesian posterior estimates) 

Analyzing the parameter estimates 6^ and 62 (see Table[7]of Appendix C) for the two regimes, 
we find that Regime 1 has stronger dependencies on the first tree, whereas Regime 2 has stronger 
dependencies on the second tree. In particular. Regime 2 exhibits stronger conditional negative 
dependencies reflected by rotated Gumbel copulas, thus creating a more asymmetric dependence 
structure. 

In order to apply our Bayesian estimation procedure, we need to distinguish both regimes to avoid 
model identification problems. For a detailed consideration of this issue we refer the reader to 



Friihwirth-Schnatter (2001) Using our observations with regard to the strength of dependence 
in the two regimes identified by the EM- Algorithm, we define Regime 1 to correspond to weaker 
dependence on the second tree and Regime 2 to correspond to stronger dependence on the second 
tree, compared by the sum of absolute values of Kendall's r corresponding to parameters 0^'*^*^*^'^ 
and Q^^'^^^^ ^ The resulting posterior probability estimates for the hidden state variable, i.e. 

P{St = l\uT) :=f;5['^'^^^, 

r=l 

for R independent MCMC samples, are plotted as dotted points in Figure [2] These Bayesian 
estimates follow those obtained from the EM algorithm closely, showing only a little bit more 
variability. 



5.1.2. Identifying crisis regimes 

Having identified parameter switches in an R-vine copula model for our data set, we will now 
try to identify switches in the overall dependence structure. Since there is empirical evidence that 



dependence structures can change in times of crisis (c.f. Longin and Solnik (1995) Ang and Bekaert 



(2002) or Garcia and Tsafack (2011)) and since tail dependencies become more important in times 



of extremal returns, we want to select two different R-vine structures using the aforementioned 
procedures. To do so, we start with a roUing window analysis, selecting and fitting R-vine models 
to a rolling window of 100 data points. To reduce model complexity, we decide to work again with 
truncated and simplified R-vines, resulting in a sufficiently fiexible and parsimonious model. The 
copulas on the first tree were chosen to be either all Gaussian, Gumbel or survival Gumbel. The 
copulas on the second tree were set to Gaussian and the R-vines were truncated after this second 
tree. The resulting rolling log likelihoods from conducting this analysis are given in Figure |3] Note 
that for AIC (BIG) comparison this is sufficient since the number of parameters remains the same 
in all models considered. 



log likelihood Normal - survival Gumbel 




200 400 600 800 200 400 600 800 

starting observation of roiling window starting observation of roiling window 

Figure 3 Left-panel: log likelihood values resulting from fitting R-vines with normal (solid), survival 
Gumbel (dashed) and Gumbel (dotted) copulas. Right panel: difference between the values for the 
normal and Gumbel model, we indicate periods to which the structures for Model (3) are fitted. 



We see that, while the range of overall likelihood estimates is similar, the Gaussian model tends 
to give the best fit, i.e. highest log likelihood, (left panel of Figure |3]) over the whole data set. 
However, the survival Gumbel model starts to outperform the Gaussian model towards the end of 
the observation period (right panel of Figure [s]). Furthermore, the survival Gumbel model, in which 
the exchange rates taken into consideration are assumed to be lower tail dependent, tends to outper- 
form the model with standard Gumbel copulas, corresponding to upper tail dependence. This is in 
accordance with the observation that the financial crisis during the observation period originated in 
the dollar area, quickly spreading to the world economy but with different severity e.g. to the other 
developed countries and the developing countries. Because of this, cash fiows out of the dollar area, 
resulting in higher FC/US exchange rates, tend to be less extremely correlated than cash fiows into 
the dollar area to settle liabilities denominating in US dollar, which results in more lower than upper 
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tail dependence. Given these observations from the rolling window analysis, we decide to select R- 
vine copulas as follows: For Regime 1, the tree structure (Vi) is again fitted to the whole data set, but 
we use only Gauss copulas as bivariate building blocks. Since parameter estimates on the higher trees 
indicate week dependence or even conditional independencies, we truncate the R-vine copula after 
the second tree. To determine a second structure (V2) we apply our outlined techniques to the time 
frame from July 10, 2008 to December 3, 2008, where the likelihood of the survival Gumbel model 
starts to be higher than the likelihood of the Gaussian model. Doing so, we also capture many high- 
impact events of the financial crisis. For the copulas associated to the R-vine structure we consider 

Model (2a) survival Gumbel copulas on the first tree to capture strong dependencies for nega- 
tive returns. 

Model (2b) Gumbel copulas on the first tree to capture dependencies in the upper tail. 

Model (2c) Student-t copulas on the first tree to cover symmetric tail dependencies. 

The copulas corresponding to edges on the second tree are again chosen to be Gaussian and we 
truncate after Tree 2. While the survival Gumbel model is preferred in the rolling window analysis, 
we include Models (2b) and (2c) to investigate the impact of different tail dependencies. Note that 
Model (2c) where we use Student-t copulas in the second structure is close to considering an R-vine 
model where the structure is kept fixed and only the parameters are subject to regime switches if 
the degrees of freedom parameters are high. In this case, we are left with two different (truncated) 
R-vine structures where all copulas are Gaussian. Since an R-vine where all bivariate building blocks 
are Gaussian leads to a multivariate Gaussian copula, the possible dependencies in the two regimes 
do only differ because of the truncation offer the second tree. Given the fact that limiting the range 
of copula families which can be associated to the R-vine as in Models (2a) - (2c) necessarily leads to 
lower AIC values we do further include a Model (3) which also has different dependence structures 
in both regimes but where all bivariate copulas are selected using the stepwise AIC criteria. To 
be more precise, the R-vine structure V3 together with the corresponding copulas for the "crisis" 
regime is selected by applying the stepwise selection procedure to the part of our dataset where 
the rotated Gumbel copula is outperforming the normal copula in the scenario analysis, which 
corresponds approximately to the last 250 observations (annotated with "crisis" in the right panel 
of Figure [3|. The R-vine structure for the "normal" regime, again with corresponding copulas is 
identified from the remaining data points, it is the same as for the "normal" regime in Models (2a) 
- (2c), Vi. While Models (2a) - (2c) are designed with specific dependence structures in order to 
investigate the changes of dependence present in the data more closely. Model (3) makes use of the 
full modeling flexibility of MS R-vines in order to provide a close fit to the data at hand. We employ 
the stepwise EM-procedure to fit MS models with the selected regimes to our data set. The resulting 
smoothed probabilities inferred from the hidden state variable for being in the non-Gaussian regime 
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using Models (2a) - (2c) are given in Figure |4] 

While the overall strength of dependence modeled in the two regimes (judging by the fitted values 
of Kendall's r, see Tables [s] and [9] of Appendix [C]) is similar for all three models, the results for 
Model (2c) with Student-t copulas are close to the results of Model (1), whereas the other two 
differ significantly. This was expected, since the model with t-copulas is close to a Gaussian copula 
model with regime switching parameters. Analyzing the estimated Kendall's r further, we find 
that in Model (2c) the Kendall's r between the Japan-US and the India-US exchange rate indicates 
negative dependence (see Table [s]). Since Gumbel and survival Gumbel copula models exhibit only 
positive dependence, this cannot be captured in Model (2a) or (2b), respectively. Replacing the 
copula for this bivariate margin by a Gauss copula (we refer to the resulting modification of model 
(2a) as (2a*)) so that it captures the negative dependence does however not significantly change the 
posterior estimates for the hidden state variable. This means that the observed difference in the 
behavior of Models (2a) and (2b) as compared to Model (2c) cannot be explained by the lack of 
Gumbel and Gumbel survival copulas to allow for negative dependence. Instead these models tend 
to be preferred during specific times of high impact events of the financial crisis, where the bivariate 
dependence structures are closer to the dependence structure of a Gumbel copula, as indicated in 
Figure [4] where some important events are annotated. 

While obtaining a sample from the posterior distribution of the latent state variable is a byproduct 
from the Bayesian estimation procedure, its main purpose is to characterize the joint posterior 
distribution of the model parameters in the R-vine copula. Figure [5] shows histogram plots of 
several marginal posterior densities for the copula parameters in the crisis regime in Model (2a*). 
As we can see, the parameter value of r/jv/j-jpy = which would correspond to independence is 
nowhere near a 90 or 95 percent CI, the dependence is significantly negative. For the copula between 
Brazil-US and China-US in contrast, the parameter values in our posterior sample are all close to 
0, which means that the two time series are only weekly dependent or maybe even independent. 
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Figure 4 Smoothed probabilities that the hidden state variable indicates the non-Gaussian regime. 
From top to bottom: model (2a) with Gumbel survival copulas, model (2b) with Gumbel copulas 
and model (2c) with Student-t copulas. The solid lines correspond to EM estimates while Bayesian 
MCMC estimates are dotted, high impact events of the financial crisis are annotated in the upmost 
graph. 
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INR - JPY JPY - CHF CHF - EUR EUR - GBP 




-0.35 -0.25 -0.15 0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.45 0.50 0.55 0.60 0.30 0.35 0.40 0.45 0.50 0.55 0.60 



CAD - GBP GBP - AUS AUS - BRL BRL - CNY 




0.25 0.30 0.35 0.40 0.45 0.50 0.25 0.30 0.35 0.40 0.45 0.50 0.55 0.35 0.40 0.45 0.50 0.55 0.60 0.00 0.05 0.10 0.15 0.20 



Figure 5 Estimated marginal posterior densities in the crisis regime of Model (2a*) with 90% CIs. 
The plotted densities correspond to the (unconditional) copulas associated to Tree 1 of the vine V2. 

Model (1) (2a) (2a*) (2b) (2c) (3) no MS 
DIG -4398 -4280 -4312 -4199 -4346 -4430 -4146 

Table 3 DIG values for the different (regime switching) R-vine models that have been considered. 
Lower values indicate a better fit of the model to the data. 



5.1.3. Model comparison 

Having discussed the stylized features of the different R-vine models we have taken into consideration 
for the exchange rate data, we want to compare them in terms of their fit to the observations at 
hand. For this, we rely on in-sample methods, and use our Bayesian Gibbs sampling procedure 



to calculate the deviance information criterion (DIG) which has been proposed by Spiegelhalter 



et al. (2002) Table |3] shows DIG values for all models under investigation, calculated using the 



Gibbs sampling procedure. For comparison purposes, we also include an R-vine model without MS, 
but where the vine tree structure has not been truncated after tree 2. The first two trees of this 
structure correspond to Structure Vi. 

Although the full R-vine model has 36 parameters and the MS R-vine models where we use truncated 
vines and one parametric pair copulas only have 32, even the worst MS R-vine copula outperforms 
the model without Markov structure, which clearly supports the use of models with time varying 
dependence in this context. The DIG values further show that in terms of in-sample fit, the model 
with standard Gumbel copulas in the crisis regimeis outperformed by the other models, which 
was to be expected from the rolling window analysis. Since the copulas in Model (1) were chosen 
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maximizing pairwise AIC, it outperforms the models where we restricted the choice of copulas. The 
best-performing model however is Model (3), where the copula families were chosen using pair-wise 
AIC but the R-vine structure differs between the regimes. This shows that MS models for all 
components of the dependence structure are more suitable for this kind of data than models where 
only the copula parameters are varying over time. 



5.2. Eurozone country indices 

To illustrate to the interested reader that MS can be detected in the dependence structure of financial 
data sets of various kind, we briefly outline two more applications. Our next data set consists of 
daily log-returns for 5 stock indices of the eurozone: the German DAX, the French CAC 40, the 
Dutch AEX, the Spanish IBEX 35 and the Italian FTSE MIB. We consider 985 observations between 
May 22, 2006 and April 29, 2010. The data has already been analyzed by Brechmann and Czado 



(2011) we refer to their Appendix A for the determination of marginal time series models. Since 



we only want to investigate the presence of MS, we consider the simplest possible model, having 
the same tree structure and copula families in both regimes. With the same technique as in Section 



5.1.1 (Model (1)) we fit this R-vine with only MS-parameters to the data (see Figure 11 Appendix 



O. Here, the two regimes differ significantly in terms of the dependency strength they describe (see 



Table 12, Appendix jCj) . While the copula parameters in the first regime correspond to values of 
Kendall's r of 0.6 — 0.75, the second regime has values between 0.75 and 0.85. For comparison, the 
probabilities for the high dependence regime are plotted in Figure [6] together with the probabilities 
for the "crisis" regime in Model (2a) for the exchange rate data and the quoted values of the STOXX 
50. It shows, that the high dependence regime for eurozone country indices becomes more relevant 
during times of a weakening economy and that it's presence is positively correlated with the presence 
of the "crisis" regime in Model (2a) of the exchange rate data, reflecting the interrelations between 
foreign exchange markets and stock markets. 



5.3. Selected German Stocks 

As a third application, we consider daily log-returns of the 10 stocks in the German stock index DAX 
with the highest market capitalization on October 25, 2010, namely: Allianz (ALV), BASF (BAS), 
Bayer (BAY), Daimler (DAI), Deutsche Bank (DBK), Deutsche Telekom (DTE), E.ON (EOA), 
RWE, SAP, Siemens (SIE). The observed time period ranges from January 17, 2001 to November 
4, 2010, resulting in 2494 daily observations. The analysis of the one dimensional margins and the 
transformation to copula data has been performed in Stober and Czado (2011)] For demonstration 



purposes, we again choose the simplest possible model. Employing again the technique of Sections 



5.1.1 and 5.2, we fit an R-vine with only MS-parameters (Figure 12, Appendix O. As for the 



indices, we identify one regime with weaker (Kendalls r on the first tree ~ 0.4) and one regime with 
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Figure 6 Smoothed probabilities (EM algorithm) that the latent variable indicates the presence 
of the high dependence regime in the STOXX data set (solid) vs. smoothed probabilities that the 
latent state variable indicates the presence of the crisis regime in the exchange rate data set (dotted) 
vs. quoted STOXX 50 values (dashed, right y-axis). 



stronger (r ~ 0.6) dependencies (see Table [Tsj Appendix [C] The probability over time that this 
second regime is present is plotted in Figure [7] It shows that the probability for the high dependence 
regime and the absolute index values are clearly negatively correlated for the period from 2001 to 
2008. Since the beginning of 2009 we observe that, despite rising index values, also the probabilities 
are in an upward trend. The negative correlation at the beginning suggests an interpretation of the 
high dependence regime as a kind of "crisis" regime, which governs the dependence during times 
of a weak economy. In this perception, the upward trend in the state probability since 2009 is an 
indicator for the remaining of the German economy in a "crisis" state and a growing destabilization 
due to the public debt problems in the eurozone. 



6. Discussion and Outlook 

This paper provides a detailed investigation of estimation methods for Markov switching regular 
vine models, which constitutes a significant contribution towards making time-varying copula mod- 
els a standard tool for the description of multivariate time series. Allowing flexible pair copula 
constructions, which can account for the stylized facts observed in the dependence among financial 
time series, to vary over time, we create a model which is limited mainly by its need for efficient 
computational treatment. The quick EM algorithm, based on the step-wise estimator for R-vine 
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Figure 7 Smoothed probabilities (EM algorithm) that the high dependence regime is present in 
the German stock data (solid, left y-axis) vs. quoted values of the DAX (dotted, right y-axis). 



copulas, which we have introduced, allows to perform ML inference in almost arbitrary dimensions. 
For a more thorough study of the resulting model parameter estimates and their uncertainty, we 
further introduced a Bayesian inference procedure which has been shown to correctly capture the 
true model in simulated data. 

Our findings in empirical applications illustrate that regime switches are present in many financial 
data sets, including exchange rates, index returns and stock returns. In particular, regime switching 
models constitute possible tools for the accurate description of changes in dependence during times 
of crisis. In this context we have also demonstrated a possible model selection heuristic for the 
context of regime switching dependence models using the R-vine selection procedures pioneered by 
Di£mann et al. (2011)| 

While we believe that the methods for parameter estimation presented here will satisfy most statis- 
ticians and practitioners needs, improvements in model selection are still desirable. It is computa- 
tionally intensive to perform model selection by computing the DIG from a sufficiently large sample 
of the parameters' posterior distribution, therefore model selection will continue to be a topic of on- 
going research. In particular, we will investigate how the rolling window methods can be improved 
and automatized using change point detection. 
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Appendices 

A. Step-wise estimation of R-vine copula parameters 

In this appendix, we will shortly describe the step wise parameter estimation method for an R-vine 
copula for the case where observations are subject to different weights. Here, the function to be 
maximized is a weighted log likelihood 

T 

^o^tW ■=Y.^og{c{nt\V,B,G))-ujt, 

t=i 

where {vit)t=i,...,T is a series of independent observations and {wt, t = 1, . . . , T} are the corresponding 
weights. 

For the case where weights are constant, this estimation procedure has already been applied by 



several authors among which we want to particularly note Haff (2010) and Czado et al. (2011) For 



more details, in particular on asymptotic properties and further references we refer to the paper 



of Haff (2010) In the context of MS models however, the weights to be attributed to individual 
observations are determined in the Expectation Step of the EM algorithm and are non-constant 
therefore, in particular ujt = {ytt\x{0^)^ ^ at iteration I. 

Before we can come to the algorithm itself, we need to introduce some more notation. Evaluating 
the density of an R-vine copula using Expression ([T]) involves arguments 

. I . _ ^Cx,x,|x_, {F{x\^-j). F{xj\yi^j)) 
aF(x,-|x_,-) 

where x_j is short notation for the vector where the j-th component of x {xj) has been removed. 
Following Aas et al. (2009), we will denote this conditional (copula) distribution function in the 
bivariate case when Xi = Ui and X2 = U2 are uniform, i.e. jointly distributed according to some 
copula C, by h{u; v , r] , t) . Then, 

h{ui;u2,V,t> = F{ui\u2) = -^-^ , 

0U2 

where the second argument U2 corresponds to the conditioning variable and (77, t) denotes the set 
of parameters rj together with the copula type t for the copula Cui,u2 of joint distribution. The 
corresponding copula density will be denoted by c{ui,U2;r],t) 

Using this, the estimation algorithm is given by the steps outlined in Algorithm [2] We use the 
convenient matrix notation for this computational purpose and our algorithm determines the correct 
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copula for each step. As in Section 2.1 we denote the parameters and copula types corresponding 



to the copula associated to edge e with 'nj{e),k{e)\D{e) ^^'^ tj(e),k{e)\D{e) respectively. 
Algorithm 2 Obtain stepwise parameter estimates for an R-vine copula. 

Require: R-vine copula specification in the form of an R-vine matrix with corresponding sets of 
copula parameters / types, and observations (ui, . . . , uj^), = (uu, ■ ■ ■ ,Udt) from the R-Vine 
copula distribution with corresponding weights {ujt, t = 1, . . . , T} 
Define mfcj = ■max{mk,i, ...,rn-„^j} 

n i. direct,t , -. ^^ ■ i j 

■=Uit t = l,...,T l = l,...,d 

for i in d-l,...,l do 
for k in d,...,i+l do 

f direct,t , -, m 

zr-^:=v^- t = l,...,T 



if m = m-fc j then 

direct,t 



f direct,t , -. 

zrn := V, , C = 1, 



else 

t indirect,t , rp 

^'^2 ■- ^fc,(d-m+l) E - 1, • • • , J 

end if 

Obtain pk^i via maximization of 

'^k-l,i • 1' ^' 2^ 'lmk,i,rni^i\mk+i^i,md^i^ '-mi:^i,mi^i\rn^.+i^i,ma^iJ -l, . . . , j. 

end for 
end for 



B. Relative bias and MSE for parameter estimates in the 
simulation study 

The following tables show the relative bias and relative MSE for the parameters of two selected 
scenarios in the simulation study (Scenario 2 & Scenario 4). For comparison purposes all copula 
parameters have been transformed to the Kendall's r level. 

Notice the large bias for the posterior mean estimate of the second Markov chain parameter in 
Scenario 4 (Table [5]) where identification issues where observed. In this case the Gibbs sampler 
with objective priors fails to capture the underlying Markov structure correctly and the Bayesian 
procedure needs to be started with strong subjective prior beliefs. 

In general, we observe that the estimation error in the second and third tree is higher than on the 
first tree and that the uncertainty in the Gumbel regime, from which less realizations are included 
in the data set, is higher than in the Gaussian regime. 
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Gaussian Regime 



— I 1 \ r~ 

43|12 42|1 32|1 41 



31 



"~r" 

21 



Gumbel Regime 



— I 1 1 — 

41|23 42|3 31 12 



21 



43 



32 



Gaussian Regime t^^i^ T4211 



"^32 1 1 



Ul 



T3I 



T2I 



relative bias 


-2.1 -10- 


-3 


-1.7 -10- 


-3 


7.2 -10- 


-4 


-8.1 -10- 


-4 


-1.1 -10- 


-3 


-1.9 -IQ- 


-3 


relative MSE 


1.4 -10- 


-3 


4.3 10- 


-4 


4.3 -10- 


-4 


4.6 -10- 


-5 


4.1 -10- 


-5 


5.5 -10- 


-5 


Gumbel Regime 


T41|23 




"^4213 




"^3112 




T43 




T32 




7-21 




relative bias 


9.5 -10- 


-2 


5.2 10- 


-3 


1.7 -10- 


-2 


-1.1 10- 


-2 


-1.6 -10- 


-2 


-1.0 -10- 


-2 


relative MSE 


1.3 -10- 


-2 


7.2 -10- 


-3 


5.6 -10- 


-3 


3.3 10- 


-3 


3.2 -10- 


-3 


3.6 -10- 


-3 



Table 4 Scenario 2: Relative error of Kendall's r estimates (top figure) and relative bias / MSE 
for the Gaussian regime (upper table) and the Gumbel regime (lower table), respectively. 



Relative error of MC parameter estimates (Scenario 2) 



MC parameters a 



relative bias -3.7 -lO"^ -1.0 -IQ-^ 
relative MSE 1.4 -lO"^ 6.5 -lO"^ 



Relative error of MC parameter estimates (Scenario 4) 



MC parameters a 



relative bias -9.7 -lO'^ -7.2 -10"^ 
relative MSE 4.0 -lO-^ 5.1 -lO-^ 



Table 5 Relative error of Markov chain parameter estimates in Scenarios 2 (left figure) and 4 (right 
figure), and relative bias / MSE for Scenario 2 (left table) and Scenario 4 (right table). 
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Gaussian Regime 

o 
8 

-s- -s- _S_ g _8_ 




Gumbel Regime 

8 -j- 

; ; ^ ^ -0- 
1 1 


^ T T ^ ^ 

o 





— I 1 \ \ 1 1 1 1 1 1 1 r 

43|12 42|1 32|1 41 31 21 41|23 42|3 31 12 43 32 21 



Gaussian Regime t^^^i2 Ti2\i T32\i "^41 T31 T21 



relative bias 


2.4 


10-1 


1.9 


•10" 


-2 


-2.6 


•10-2 


-5.1 


•10-2 


-3.7 


•10 


-2 


-1.5 


•10^ 


-2 


relative MSE 


2.5 


10-2 


7.5 


•10" 


-3 


8.3 


10-3 


5.7 


10-3 


5.4 


•10- 


-3 


4.4 


•10- 


-3 


Gumbel Regime 


'^41123 


"^4213 




"^3112 




T32 








relative bias 


3.1 


•lOO 


8.7 


•10- 


-1 


1.2 


•10° 


5.2 


•lOi 


5.1 


•10- 


-1 


4.5 


•10- 


-1 


relative MSE 


2.3 


•10° 


5.3 


•10- 


'1 


7.5 


10-1 


2.7 


10-1 


2.3 


•10- 


-1 


2.1 


•10- 


-1 



Table 6 Scenario 4: Relative error of Kendall's r estimates (top figure) and relative bias / MSE 
for the Gaussian regime (upper table) and the Gumbel regime (lower table), respectively. 

C. Selected R-vine structures 

C.l. Selected R-vine structure for the US exchange rates 
Model (1) 



In Section 5.1.1| an R-vine is fitted to the exchange rate dataset using the outlined procedures of 



Difimann et al. (2011)| and Brechmann et al. (2011) Here only the copula parameters are switching 



while tree structure and copula families are common to both regimes. The resulting tree structure 
is given in Figure |8] the following together with corresponding Kendall's r estimates in Table [7] 
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Tree 1 


GBP,EUR EUR,CHF CHF,JPY AUS,EUR AUS,BRL INR,AUS CAD,AUS CNY,INR 


cop. fam. 


SG 


SG 


SG 


N 


G 


N 


N 


G 








Re^ 


^ime 1 










-EM 
^1 

^MCMC 

5% quant. 
95% quant. 


0.55 
0.56 
0.53 
0.60 


0.78 
0.79 
0.77 
0.81 


0.46 
0.47 
0.43 
0.50 


0.46 
0.46 
0.43 
0.49 


0.19 
0.18 
0.14 
0.23 


0.14 
0.14 
0.10 
0.20 


0.29 
0.28 
0.24 
0.32 


0.11 
0.13 
0.09 
0.17 










;ime 2 










-EM 
^2 

■^MCMC 

5% quant. 
95% quant. 


0.44 
0.44 
0.40 
0.47 


0.58 
0.58 
0.55 
0.60 


0.24 
0.22 
0.17 
0.27 


0.41 
0.40 
0.36 
0.43 


0.45 
0.43 
0.39 
0.47 


0.26 
0.25 
0.20 
0.29 


0.44 
0.44 
0.40 
0.50 


0.07 
0.05 
0.02 
0.10 



^ ^ GBP,AUS| CAD,EUR| CAD,BRL| BRL,INR| CNY,AUS| JPY,EUR| CHF,AUS| 





EUR 


AUS 


AUS 


AUS 


INR 


CHF 


EUR 


cop. fam. 


G 


G 


SG 





SG 


G270 


G270 








Regime 


1 










0.15 


0.11 


0.07 


0.02 


0.01 


-0.06 


-0.03 


jt.MCMC 
^1 


0.15 


0.11 


0.07 


0.01 


0.02 


-0.06 


-0.03 


5% quant. 


0.10 


0.07 


0.03 


-0.04 


0.00 


-0.10 


-0.07 


95% quant. 


0.20 


0.15 


0.13 


0.06 


0.05 


-0.02 


-0.00 








Regime 


2 








-EM 


0.15 


0.11 


0.11 


0.11 


0.10 


-0.31 


-0.24 


-MCMC 


0.16 


0.13 


0.11 


0.11 


0.11 


-0.31 


-0.24 


5% quant. 


0.10 


0.08 


0.04 


0.06 


0.07 


-0.36 


-0.28 


95% quant. 


0.21 


0.17 


0.17 


0.16 


0.16 


-0.26 


-0.19 



Table 7 Estimated Kendall's r for the first and second tree of Model (1). 
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jpy 




EUR.CHF GS70 EUR.AUS 



Figure 8 The first and second tree of tlie R-vine tree structure Vi of Model (1) for the exchange 
rate data. We also choose this structure for the non-crisis regime in Models (2a) - (2c) and Model 
(3). 

Model (2) 

For Models (2a) - (2c), two R-vine tree structures have been selected to account for different de- 
pendencies during times of crisis and normal times. The first R-vine (Vi), corresponding to normal 
times, has again the structure displayed in Figure [8] the first and second tree of the second R-vine 
(V2), corresponding to crisis times, is given in Figure^ 



/ 




\ 



Figure 9 The first and second tree of the "crisis" R-vine structure V2 which we have chosen for 
Model (2) of the US exchange rate data. 
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"normal", Vi GBP,EUR EUR,CHF CHF,JPY AUS,EUR AUS,BRL INR,AUS CAD,AUS CNY,INR 



(2a) f f ^ 


0.53 


0.75 


0.45 


0.46 


0.28 


0.21 


0.35 


0.12 


\-^°') ' 1 


0.56 


0.78 


0.45 


0.48 


0.24 


0.21 


0.33 


0.15 


5% quant. 


0.52 


0.75 


0.41 


0.45 


0.19 


0.17 


0.29 


0.10 


95% quant. 


0.60 


0.80 


0.49 


0.51 


0.29 


0.25 


0.36 


0.20 


(2h) ff 


0.54 


0.75 


0.44 


0.47 


0.29 


0.22 


0.35 


0.12 


(2h) ff ^""'^ 


0.52 


0.74 


0.43 


0.44 


0.29 


0.21 


0.34 


0.11 


O/O C[U.clIlX. 




u. / Z 






U.ZO 


n 1 7 

U. 1 1 


U.Ol 


n n7 

U.U 1 


95% quant. 


0.55 


0.76 


0.46 


0.48 


0.32 


0.24 


0.37 


0.14 


(2c) ff ^ 


0.60 


0.81 


0.47 


0.48 


0.21 


0.19 


0.32 


0.16 


(2c) ff 


0.60 


0.80 


0.47 


0.49 


0.21 


0.21 


0.31 


0.17 


5% quant. 


0.57 


0.79 


0.44 


0.46 


0.16 


0.17 


0.27 


0.13 


95% quant. 


0.62 


0.82 


0.51 


0.52 


0.26 


0.22 


0.35 


0.22 



"crisis", V2 GBP,EUR EUR,CHF CHF,JPY JPY,INR AUS,GBP BRL,AUS BRL,CNY CAD,GBP 



(2a) ff ^ 


0.44 


0.45 


0.11 


0.00 


0.41 


0.49 


0.11 


0.41 


(2a) ff 


0.42 


0.52 


0.22 


0.01 


0.37 


0.47 


0.07 


0.37 


5% quant. 


0.36 


0.47 


0.11 


0.00 


0.30 


0.41 


0.01 


0.29 


95% quant. 


0.49 


0.56 


0.30 


0.02 


0.44 


0.53 


0.13 


0.43 


/r,i ^ - EM 
(2b) T2 


0.37 


0.37 


0.10 


0.00 


0.32 


0.41 


0.08 


0.36 


(2b) ff 


0.45 


0.37 


0.05 


0.01 


0.35 


0.44 


0.13 


0.38 


5% quant. 


0.34 


0.27 


0.00 


0.00 


0.25 


0.36 


0.03 


0.30 


95% quant. 


0.55 


0.45 


0.15 


0.04 


0.44 


0.53 


0.24 


0.46 


(2c) f r 


0.43, 10.8 


0.58, 8.6 


0.27, 7.9 


-0.13, 30 


0.37, 10.7 


0.44, 5.9 


0.06, 30 


0.33, 30 


(2c) ff 


0.42, 14.2 


0.56, 9.7 


0.25, 9.8 


-0.15, 21.4 


0.35, 15.4 


0.45, 9.3 


0.05, 20.8 


0.34, 21.8 


5% quant. 


0.38, 7.0 


0.53, 5.6 


0.21, 5.2 


-0.21, 11.6 


0.31, 7.0 


0.45, 4.8 


-0.01, 10.7 


0.29, 11.0 


95 % quant. 


0.46, 25.7 


0.60, 16.3 


0.29, 18.3 


-0.11, 29.2 


0.40, 28.1 


0.49, 19.0 


0.11, 29.2 


0.39, 29.3 



Table 8 Estimated Kendall's r values corresponding to the first tree of Models (2a) - (2c), respec- 
tively. For the t-copula used in Model (2c), the first parameter is transformed to Kendall's r, the 
second parameter gives the estimated degrees of freedom (with = 30 as upper limit. 
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0.12 
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0.04 


0.30 
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5 % quantile 
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0.25 


0.12 


-0.20 


-0.40 


0.05 
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Table 9 Estimated Kendall's r, corresponding to the second tree of Models (2a) - (2c). 



Model (3) 

Also for Model (3), we consider two different R-vine structures too account for changes in depen- 
dence. The structure for the first regime is again chosen to be Vi, with copulas selected by AIC, 
the tree structure for the second regime is V3, selected from the part of the data labeled "crisis" in 
Figure |3] and shown in Figure 10 The corresponding copulas, together with estimated values for 



Kendall's r, are given in Table 10 for the "normal" regime and in Table 11 for the "crisis" regime, 
respectively. 
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AUS.BRL 

EUR.AUS 



Figure 10 The first and second tree of tlie R-vine structure which we have chosen for the "crisis" 



regime of Model (3) . 


We refer to this structure as V3. 
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Table 10 Estimated Kendall's r, corresponding to the first and second tree of the normal regime 
in Model (3). 
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0.19 


0.22 


0.41 


-0.27 



Table 11 Estimated Kendall's r, corresponding to the first and second tree of the crisis regime in 
Model (3). 



C.2. Selected R-vine structure for Eurozone country indices 

For the return data set of Eurozone country indices we selected a common R-vine structure (Figure 



11), with common copula families (Table 12) but MS parameters. The estimated Kendall's r values 



are given in Table 12 




CACaO.FTSE MIB 



Figure 11 First and second tree of the R-vine structure selected for Eurozone country indices. 
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Tree 1 



DAX,CAC40 IBEX,CAC40 AEX,CAC40 FTSE MIB,CAC40 



cop. fam. 


SG 


N 


N 


SG 


Ti (Kegime 1) 


0.71 


0.62 


0.72 


0.66 


T2 (Hegime 2) 


0.83 


0.78 


0.81 


0.77 



AEX,DAX I DAX,FTSE MIB I FTSE MIB,IBEX I 
Tree 2 > 1 > 1 ' 1 

CAC40 CAC40 CAC40 

cop. fam. N G N 

Tf^(Regimel) 0.06 0.13 0.12 

rf*^ (Regime 2) 0.17 0.15 0.17 

Table 12 Values of Kendall's r estimated for the county indices. 
C.3. Selected R-vine structure for the German stock index (DAX) 

Treel 




Figure 12 First and second tree of the R-vine structure selected for German stock returns. The 
corresponding values of Kendall's r are given in Table 



13 



Tree 1 DTE,ALV ALV,DBK DBK,SIE SIE,DAI SIE,SAP SIE,BAS BAS,BAY BAS,EOA EOA,RWE 



cop. fam. 


SG 


SG 


SG 


SG 




SG 


SG 


SG 


SG 


G 


Ti (Kegime 1) 


0.26 


0.42 


0.38 


0.34 




0.34 


0.32 


0.39 


0.25 


0.49 


T2 (Kegime 2) 


0.54 


0.63 


0.64 


0.63 




0.59 


0.65 


0.61 


0.55 


0.66 


Tree 2 


DTE,DBK| 
ALV 


SIE,ALV| SAP,DBK 
DBK SIE 


DAI,DBK| DAI,BAS| 
SIE SIE 


SIE,BAY1 
BAS 


E0A,SIE1 
BAS 


BAS,RWE 
EOA 


copula 


G 


G 


G 




G 




G 


G 


G 


N 


Ti (Kegime 1} 


0.16 


0.16 


0.19 




0.22 




0.19 


0.10 


0.12 


0.07 


-EMf^ . 

T2 (Kegime 2) 


0.22 


0.33 


0.13 




0.23 




0.20 


0.25 


0.25 


0.18 



Table 13 Values of Kendall's r determined for the DAX data set using the EM algorithm. 



36 



References 



Aas, K., C. Czado, A. Frigessi, and H. Bakkcn (2009). Pair-copula construction of multiple 
dependence. Insurance: Mathematics and Economics 44i 182-198. 

Albert, J. H. and S. Chib (1993). Bayes inference via Gibbs sampling of autoregressive time series 
subject to Markov mean and variance shifts. Journal of Business & Economic Statistics 
pp. 1-15. 

Almeida, C. and C. Czado (2011). Efficient bayesian inference for stochastic time- varying copula 
models. Computational Statistics and Data Analysis. To appear. 

Ang, A. and G. Bekaert (2002). International asset allocation with regime switching. Review of 
Financial Studies 11, 1137-1187. 

Ang, A. and J. Chen (2002). Asymmetric correlations of equity portfolios. Journal of Financial 
Economics 63, 443-494. 

Bedford, T. and R. Cooke (2001). Probability density decomposition for conditionally dependent 
random variables modeled by vines. Ann. Math. Artif. Intell. 32, 245-268. 

Bedford, T. and R. Cooke (2002). Vines - a new graphical model for dependent random variables. 
Annals of Statistics 30, 1031-1068. 

Brechmann, E. and C. Czado (2011). Extending the CAPM using pair copulas: The regular vine 
market sector model, preprint. 

Brechmann, E., C. Czado, and K. Aas (2011). Truncated regular vines in high dimensions with 
applications to financial data, preprint. 

Carlin, B. P. and T. A. Louis (2009). Bayesian Methods for Data Analysis (third ed.). Chapman 
L Hall/CRC, Boca Raton, Florida. 

ChoUete, L., A. Heinen, and A. Valdesogo (2009). Modeling international financial returns with 
a multivariate regime-switching copula. Journal of Financial Econometrics 7(4), 437-480. 

Czado, C. (2010). Pair-Copula Constructions of Multivariate Copulas. In P. Jaworski, F. Durante, 
W. K. Hardle, and T. Rychlik (Eds.), Copula Theory and Its Applications, Lecture Notes in 
Statistics, Volume 198, Heidelberg, pp. 93-109. Springer- Verlag. 

Czado, C, F. Gartner, and A. Min (2010). Joint Bayesian inference of D-vines with AR(1) 
margins. In H. Kurowicka, D. Joe (Ed.), Dependence Modeling-Handbook on Vine Copulas, 
Singapore, pp. 330-359. World Scientific Publishing. 

Czado, C, U. Schepsmeier, and A. Min (2011). Maximum likelihood estimation of mixed C-vines 
with application to exchange rates. Statistical Modelling. To appear. 



37 



Difimann, J. (2010). Statistical Inference for Regular Vines and Application. Diplomarbeit, Center 
of Mathematical Sciences, Munich University of Technology, Garching bei Miinchen. Available 
under http:/ /www-m4.ma.tum.de/Diplarb/da_txt.html. 

Difimann, J., E. C. Brechmann, C. Czado, and D. Kurowicka (2011). Selecting and estimating 
regular vine copulae and application to financial returns, submitted. 

Engle, R. F. (1982). Autoregressive conditional heteroscedasticity with estimates of the variance 
of united kingdom inflation. Econometrica 50(4), 987-1007. 

Filardo, A. J. and S. F. Gordon (1998). Business cycle durations. Journal of Econometrics 85 {1), 
99.123. 

Friihwirth-Schnatter, S. (2001). Markov Chain Monte Carlo estimation of classical and dynamic 
switching and mixture models. Journal of the American Statistical Association 5^(453), 194- 
209. 

Gamerman, D. and H. F. Lopes (2006). Markov Chain Monte Carlo - Stochastic Simulation for 
Bayesian Inference. Chapman &; Hall / CRC. 

Garcia, R. and G. Tsafack (2011, August). Dependence structure and extreme comovements in 
international equity and bond markets. Journal of Banking & Finance 55(8), 1954-1970. 

Haff, I. H. (2010). Estimating the parameters of a pair copula construction, preprint. 

Hafner, C. M. and H. Manner (2010). Dynamic stochastic copula models: estimation, inference 
and applications. Journal of Applied Econometrics. 

Hamilton, J. D. (1989). A new approach to the economic analysis of nonstationary time series 
and the business cycle. Econometrica 57, 357-384. 

Hamilton, J. D. (1990). Analysis of time series subject to regime changes. Journal of Economet- 
rics 45, 39-70. 

Heinen, A. and A. Valdesogo (2009a, November). Asymmetric capm dependence for large dimen- 
sions: the canonical vine autoregressive model. CORE Discussion Papers 2009069, Universite 
catholique de Louvain, Center for Operations Research and Econometrics (CORE). 

Heinen, A. and A. Valdesogo (2009b). Truncating canonical vine copulas. Research Paper. 

Hofmann, M. and C. Czado (2010). Assesing the VaR of a portfolio using D-vine copula based 
multivariate GARCH models, in revision. 

Joe, H. (1996). Families of m-variate distributions with given margins and m(m-l)/2 bivariate 
dependence parameters. In L. Riischendorf and B. Schweizer and M. D. Taylor (Ed.), Dis- 
tributions with Fixed Marginals and Related Topics, Volume 28, Hayward, CA, pp. 120-141. 
Inst. Math. Statist. 



38 



Joe, H., H. Li, and A. K. Nikoloulopoulos (2010, January). Tail dependence functions and vine 
copulas. Journal of Multivariate Analysis 101(1), 252-270. 

Joe, H. and J. J. Xu (1996). The estimation method of inference functions for margins for mul- 
tivariate models. UBC, Dept. of Statistics, Technical Report 166. 

Kass, R. E., B. P. Carlin, A. Gelman, and R. M. Neal (1998). Markov Chain Monte Carlo in 
practice: A roundtable discussion. The American Statistician 52 (2), 93-100. 

Kim, C.-J. and C. R. Nelson (1998). Business cycle turning points, a new coincident index, and 
tests of duration dependence based on a dynamic factor model with regime switching. The 
Review of Economics and Statistics 80(2), 188-201. 

Kotz, S., N. Balakrishnan, and N. L. Johnson (2000). Continuous Multivariate Distributions, 
Volume 1: Models and Applications (Second ed.). John Wiley &; Sons, New York. 

Lewandowski, D., D. Kurowicka, and H. Joe (2009, October). Generating random correlation 
matrices based on vines and extended onion method. J. Multivar. Anal. 100, 1989-2001. 

Longin, F. and B. Solnik (1995). Is the correlation in international equity returns constant: 1960- 
1990? Journal of International Money and Finance -Z-^(l), 3-26. 

Longin, F. and B. Solnik (2001). Extreme correlations in international equity markets. Journal 
of Finance 56, 649-676. 

Manner, H. and O. Reznikova (2011). A survey on time- varying copulas: Specification, simulations 
and estimation, forthcoming in Econometric Reviews. 

Markwat, T., E. Kole, and D. van Dijk (2009, November). Contagion as a domino efi^ect in global 
stock markets. Journal of Banking & Finance 55(11), 1996-2012. 

Min, A. and C. Czado (2010, Fall). Bayesian inference for multivariate copulas using pair-copula 
constructions. Journal of Financial Econometrics 8(4), 511-546. 

Morales-Napoles, O. (2008). Bayesian Belief Nets and Vines in aviation safety and other appli- 
cations. Ph. D. thesis, Technische Universiteit Delft. 

Okimoto, T. (2008). New evidence of asymmetric dependence structures in international equity 
markets. Journal of Financial and Quantitative Analysis ^5(03), 787-815. 

Pelletier, D. (2006). Regime switching for dynamic correlations. Journal of Econometrics 131(1- 
2), 445-473. 

Prim, R. C. (1957). Shortest connection networks and some generalizations. Bell System Tech- 
nology Journal 36, 1389-1401. 

Spiegelhalter, D. J., N. G. Best, B. P. Carlin, and A. Van Der Linde (2002). Bayesian measures 
of model complexity and fit. Journal of the Royal Statistical Society: Series B (Statistical 
Methodology) 64.(A), 583-639. 



39 



Stober, J. and C. Czado (2011). Sampling pair copula constructions with applications to mathe- 
matical finance, preprint. 

Valdesogo, A. (2009). Multivariate volatility models using copulas. Ph. D. thesis, Center for Op- 
erations Research and Econometrics, Universite Catholique du Louvain. 

Wu, C. F. J. (1983). On the convergence properties of the em algorithm. The Annals of Statis- 
tics pp. 95-103. 



40 



